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Abstract 

We analyze the I-V characteristic of a double quantum well integrated in 
an antenna and driven by THz radiation. We propose a microscopic model 
to calculate the sequential current including selfconsistently the Coulomb in- 
teraction in a mean held approximation. Additional features in the current 
appear which depend not only on the frequency and intensity of the ac held 
but also on the charge in the wells. We observe, for some intensities and 
frequencies of the ac held, dynamical localization, i.e., quenching of the trans- 
mission, and, in the linear response regime, absolute negative conductance. 
We compare our calculations with recent experiments in driven double wells. 
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The effect of an ac field on the transport properties of nanostructures 
has been a subject of increasing interest in the last years. The experiments 
on tunneling through a metal-superconductor diode [1] irradiated with mi- 
crowaves indicated that the external field assists the tunneling openning new 
effective channels. Recently transport experiments have been performed in 
semiconductor nanostructures in different ac fields: The photoassisted cur- 
rent through a double barrier (DB) irradiated with FIR light shows [2,3] 
additional features which can be explained in terms of the coupling of differ- 
ent electronic states induced by the electromagnetic field [4]. The effect of an 
AC potential on the conductance and current in quantum dots (QD's) [5-7] 
, superlattices (SL's) [8] and double quantum wells (DQW's) [9] has been 
measured and new physical features as dynamical localization and absolute 
negative conductance have been observed [8]. Also Rabi transitions between 
discrete states and electron pumping in double quantum dots (DQD's) has 
been analyzed within the Keldysh formalism [10]. In this work we propose a 
model for analyzing the time-average current in a triple barrier integrated in 
an antenna and driven by THz radiation. In this configuration, the irradiated 
antenna produces an oscillatory signal between the left and right lead, i.e., a 
time dependent bias drops between the emitter and collector. It is important 
to stress that this configuration of the ac field is different than the one corre- 
sponding to homogeneous irradiation of the sample [2-4]. In order to account 
with the scattering due to roughness at the interfaces, phonons or impurities 
we analyze sequential tunnel: we calculate the current from the emitter to 
the left well J e< i, from this one to the right well an d from the right well 
to the collector J2,c- Current conservation determinates the Fermi energies 
in the wells and the sequential current. The effect of electron-electron inter- 
action in these quasi-3D structures can be treated as a perturbation to the 
non-interacting system within a mean field approximation [11,12]. The charge 
accumulated in the wells produces an additional electrostatic potential which 
modifies their electronic structure, transmission coefficient and current. We 
solve selfconsistently the Schrodinger equation including the ac field and the 
Poisson equation. The induced electrostatic potential does not drop linearly 
and accumulation and depletion layers in the emitter and collector are built 
up. Details of the electrostatic model are given elsewhere [12]. 
The hamiltonian of the DQW driven by an ac potential is: 
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Here T^kj , Tfcifc 2 are ^ ne tunneling matrix elements (leads- well and well- well 
respectively) [13]. The ac field induces a time dependent bias between the 
emitter and collector producing a dipole around the central region. The elec- 
tronic energies become time dependent Ek t {t) — E®. + eFziCosu>ot, where 
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Zi is taken constant (its mean position) in each spatial region. It implies 
that the mixing of electronic states within each spatial region due to the 
position operator Z is neglected. This approximation, which was assumed 
in the model of Tien and Gordon [14] is reasonable for this configuration 
of the ac field, where the time dependent potential drops just between the 
leads and of course in non resonant conditions, i.e, for loq ^ ^ (A be- 
ing the subband spacing) . The operators acquire phase factors of the form 
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In the case of sequential tunneling first order in perturbation theory gives a 
good description of the tunneling (the coherent tunneling, i.e, infinite order 
in the tunneling perturbation has been solved in [13]). We obtain for the 
transition probabilities: 
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Where E^ — ^— If + e^ z , the rest of parameters are shown in fig.l, /? e i — 

eF( ~%LoJ^ ' — eF( ~%LoJ^ ani ^ ^ 2c — eF ^n^ Z2 ^ are reia ted with the ac potential 
drops between the different regions, and the terms in brackets are the inelastic 
transmissions Ti(e e , ei), T 2 (t\, e 2 ), T 3 (e 2 , e c) through each barrier. The time 
describing the relaxation between the states due to coupling with phonons 
is ~ 10 _13 s [15]. The DOS in the growth direction is described then as 
Lorentzians of 1 meV half width. In formula (2) irreversibility is assumed. In 
this case the Fermi Golden Rule is valid [16]: it is clearly the case for Pi and P3 
where there is a continuum of states involved. In the case of the transmission 
probability from the left to right well we assume irreversibility due to the 
broadening of the DOS induced by scattering as well as to the 2D continuous 
in-plane DOS. On top on that, if the tunneling time associated to the outer 
barrier is shorter than the one corresponding to the barrier separating the 
two wells the process of tunnel becomes irreversible. However, in the case of 
an isolated DQD, non connected with a continuous DOS at the leads and in 
the ideal case were no scattering occurs the tunnel is a reversible process and 
a description in terms of Rabi transitions should be considered [10,16]. The 
expression for the currents is: 
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Where e r \ and e r2 are the ground states of the left and right wells. For sim- 
plicity we restrict the previous equations to the tunneling between ground 
to ground state (the generalization to excited well states is straightforward). 
The Fermi energies e w \ and t w i of the wells are obtained through the set of 
transcendental equations J — J e \ — 0; J — J\2 — 0; J — J 2 c = 0. At this point 
we solve the Poisson and Schrodinger equation selfconsistently to obtain the 
current. In order to compare with recent experiments [9], we analyze a sam- 
ple consisting in two GaAs quantum wells 180 and 100A wide separated by 
a barrier of Ga^Al^As of 60A. The outer barriers are 35A thick (sample a). 
The emitter and collector are n + doped (10 18 cm -3 ). In figure (2. a) the J/V 
curve is shown with and without Coulomb interaction. We observe a main 
peak coming from the alignement of the two ground states of the two wells. 
Two additional peaks show up at both sides of the main peak symmetric in 
bias but asymmetric in intensity and correspond to induced one photon ab- 
sorption and emission. The asymmetry in intensity takes place due to the 
increase in the transmission probability at higher bias and to the asymmetry 
of the charge accumulated in the wells (see fig. 2.b). The behaviour of J 
for different F and fixed loq is shown in fig. 3. a: as F increases new satel- 
lites show up indicating that multiphoton processes increase in probability, 
the main peak is reduced and the satellites become more intense. Also J is 
quenched for some frequencies and intensities of the field. This effect has been 
discussed [13,17,18] in terms of zeros of the Jo Bessel functions: as the ar- 
gument approaches the first zero ( ~ 2.4) dynamical localization of electrons 
takes place, independently of the transparence of the barriers. The exper- 
iments [9] are very similar to our results. However we did not observe the 
small shift that the peaks experience as the intensity increases and which is 
interpretated as heating of the contact regions [9]. As the power increases, 
the effect of the ac field on the charge accumulated in the wells increases too. 
It would be expected [11] that for high intensities the ac field would affect the 
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charge and the selfconsistent current explaining also the shift of the peaks. 
We do not observe an appreciable shift for intensities of the order of 10 5 ^- but 
it is observed for intensities at least one order of magnitude higher (as in the 
DB case [11]). There is also another effect which could explain this shift con- 
sidering the dipole term as a quantum operator [4,19]. The term eFZicos(uiot) 
will produce a shift of the energies: 8 — e F ^j Z ^ • Due to the smallness 
of the diagonal matrix elements of Z the shift becomes negligible [4]. J/V for 
different frequencies is shown in fig. 3.b. As loq increases the satellite peaks 
move far appart linearly from the main one which increases in intensity as the 
satellite peaks decrease. Finally we observe absolute negative conductance in 
a different sample: a Ga^Al^As triple barrier of 50A and a DQW of GaAs of 
150A and n + — 10 17 cm -3 (sample b). The reason for that is the asymmetry 
induced by the harmonic term Fzcosu>ot in the hamiltonian. This effect which 
has been discussed for different systems [17,20] has been observed in SL's [8] 
in the low bias region. The negative conductance induced by the ac field can 
be explained in terms of the Fermi distributions. Equation (3) shows that the 
supply function contains a superposition of effective dc voltages ^P^ each 
one weighted nonlinearly by the Bessel functions. This effect is responsible of 
the negative conductance because the Fermi factors can be expressed in terms 
of effective Fermi energies (for instance, Ep effec — Ep — eV ± pfowo, in the 
collector). For some ranges of F and loq the electrons are able to overcome 
the static bias (eV < p ma x^o) an d electronic pumping at zero bias occurs 
[21]. We study the dependence of the absolute negative conductance with F 
including also the effect of Coulomb interaction (see fig. 4a, b). We observe 
that the charge interaction changes in some cases the sign of the differential 
conductance and the current. Therefore it modifies qualitatively the current 
in the low bias regime and has to be included to describe it properly. Also 
we observe that tunning F the differential conductance changes sign as well 
(fig. 4a) as the experimental evidence shows. Also in agreement with the 
experiments the conductance becomes negative at low bias increasing F and 
again positive as F increases further [8]. 

In conclusion, the ac sequential current through a triple barrier is evalu- 
ated including the Coulomb interaction, scattering effects and temperature. 
By means of our model we can explain the behaviour of the current versus 
voltage as a function of the ac frequency and intensity in good agreement with 
the experimental information. We analyze different samples and we observe 
interesting features in the current as dynamical localization and absolute neg- 
ative conductance. Both effects can be controlled by tailoring the sample and 
harmonic potential configuration. 
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FIGURES 



FIG. 1. Schematic diagram of sequential tunneling through a triple barrier, in an ac potential. 

FIG. 2. a) J/V for ac sequential tunneling through a triple barrier (with and with- 
out selfconsistency) (inset without ac); b) 2D selfconsistent charge in the left and right 
well as a function of V (the inset shows the Fermi levels of the two wells ) for sample 
a.(F= 1.510 5 ^,u;o= 15THz,T= 100°K). 

FIG. 3. a) Selfconsistent J/V for an ac u>o=1.5 THz and different F ; b) at F — 1.510 5 ^ and 
different ac frequencies (sample a) (T — 100° K). 

FIG. 4. a) Selfconsistent J/V for an ac u>o=1.5 THz and different F in the low bias regime 
(sample b); b) Same as a) but without Coulomb interaction (T — 100° K). 
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